More about excited bottomonium radiative decays 
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Radiative decays of bottomonium are revisited, focusing on contributions from higher-order rela- 
tivistic effects. The leading relativistic correction to the magnetic spin-flip operator at the photon 
vertex is found to be particularly important. The combination of 0(v 4 ) effects in the nonrelativistic 
QCD action and in the transition operator moves previous lattice results for excited T decays into 
agreement with experiment. 

PACS numbers: 12.38.Gc, 13.20.Gd, 14.40.Pq 



In a recent paper [l[ we studied bottomonium radia- 
tive decays using lattice nonrelativistic QCD (NRQCD) 
0. It was shown that robust signals for transitions 
among the ground state and the first two excited states 
could be obtained using multiexponential fitting tech- 
niques [3j- The qualitative features of the decay ampli- 
tudes were in agreement with phenomenological expecta- 
tions but quantitatively the values obtained were consid- 
erably larger than those determined by experiment. 

In Ref. pj the b quarks were described using an 0(v 4 ) 
lattice NRQCD action and only the simplest magnetic 
spin-flip transition operator was used. In this work 
higher-order relativistic effects in the action and addi- 
tional terms in the transition operator are investigated. 
It is expected that the hindered transitions, which in- 
volve states with different principal quantum numbers, 
are very sensitive to relativistic corrections. In compar- 
ison to Ref. the additional effects considered in this 
work substantially reduce the amplitudes for excited T 
decays and bring them into line with experimental values. 

The general setup for this study is the same as in 
Ref. The gauge fields come from a 2+1-flavor dynam- 
ical simulation done on a 32 3 x 64 lattice by the PACS-CS 
collaboration An ensemble of 192 configurations is 
used. The light quark parameters are such that the pion 
mass is near physical at 156 MeV. Landau link tadpole 
improvement is implemented with a value of 0.8463 for 
the link. The 6-quark bare mass is 1.945 (in lattice units) 
and a stability parameter 11 of 4, in line with Ref. [5], is 
used. 

The three-point functions for vector to pseudoscalar 
transitions are constructed using a sequential source 
method. Starting with a vector (pseudoscalar) opera- 
tor at the source t s , the quark propagator is evolved to a 
time T at which a pseudoscalar (vector) operator is ap- 
plied. This quantity is then evolved backward in time. 
At intermediate times t s < t' < T a transition opera- 
tor is inserted and evolution is continued to complete the 
quark antiquark loop at the source. 

With a vector operator at the source and pseudoscalar 
at the sink the three-point function is expected to have 



the form 
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where the subscripts o, o' indicate the type of operator 
(local or smeared) that is used. As in Ref. [1], three types 
of operators are used: a local operator, a wavefunction 
smeared operator [|| and an operator with the smearing 
applied twice [f|. See Ref. pj for a description of the 
smearing function and parameters. 

The overlap coefficients c and simulation energies E 
are the same ones that appear in the two-point function. 
The quantity A^f' is the matrix element of the transi- 
tion operator between the vector state n and the pseu- 
doscalar state n' . The three-point function with pseu- 
doscalar source and vector sink has the same form with 
V and P labels reversed. The matrix elements A^J^ are 
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related to those appearing in ([T]) by A, 
matrix elements can be determined by fitting the t' de- 
pendence of the three-point function for a fixed T using 
overlap coefficients and energies obtained from a fit to 
two-point correlators. 

The calculations in Ref. [lj were done using an NRQCD 
action including terms to 0(v A ). In this work 0(v 6 ) 
terms are also considered. The complete action can be 
found, for example, in Ref. [7ij. Here we display only 
the spin-dependent terms linear in chromoelectric and 
chromomagnetic fields which are relevant for subsequent 
discussion (see Ref. 0] for detailed explanation of nota- 
tion): 
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where SH^ gives the 0(v 4 ) chromomagnetic coupling 
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Table I: Transition operators used in calculating three-point 
functions. The momentum k is chosen to have a component 
only in the one direction and Aifc 



(Aie ikx + e ikx Ai). 





magnetic 


electric 


0(v 4 ) 
0(v 6 ) 
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and spin-orbit terms and SH^ their leading relativistic 
corrections. The fields are tadpole-improved and in our 
calculations the action coefficients take the tree-level val- 
ues. The transition operator should be constructed in 
a way which is consistent with the action. This can be 
achieved by replacing in ^ and (O the SU(3) color elec- 
tric and magnetic fields by external electromagnetic fields 
and the strongcoupling constant by the charge. See for 
example Ref. 

For an Ml decay, the photon momentum, polarization 
vector and quark spin operator should be mutually or- 
thogonal. The explicit choice made for the operators 
used in our lattice simulation is given in Table [TJ The 
normalization is such that in the infinite mass limit the 
transition matrix element goes to f, matching the non- 
relativistic wavefunction overlap (see Eq. (1) in Ref. 

Results for the three-point transition matrix elements 
are given in Table |TT] for different source-sink time sep- 
arations and different values of recoil momentum. The 
calculations are done with the 0(v 4 ) action and the lead- 
ing magnetic operator 03. The values for T — t s equal to 
19 and 27 are from Ref. [lj. For the excited to ground 
state transitions, there is very good agreement between 
results using different time separations. However, using 
a smaller time separation T — t s = 15 allows for a better 
determination of the 2—^2 transition and we use only 
this for the present study. 

Changing the action from 0(v 4 ) to 0(v 6 ) results in a 
decrease in the strength of the spin-dependent interac- 
tions. This is evidenced by a decrease in the mass dif- 
ference between T and r\h states [7|. For our calculation, 
done at a single lattice spacing of about 0.09 fm the T-77& 
mass difference is 56(1) MeV and 24(3) MeV for IS and 
2S states respectively using the 0(v 4 ) action. These val- 
ues are reduced to 43(2) MeV and 14(2) MeV when 0(v 6 ) 
effects are included. It is expected that inclusion of radia- 
tive corrections to the coefficients of the NRQCD action 
would raise these values somewhat §. For reference, the 
experimental values are 69.3±2.8 MeV (from the Particle 
Data Group [nj) and 59.3±1.9±f"t MeV (from Adachi et 
al. [Illfor IS. For 2S there are claims of 48.7±2.3±2.1 
MeVllJl and 24.3±\ . 5 MeV [3 for the spin splitting. 

Decreasing the strength of the spin-dependent inter- 
actions leads to greater overlap of wavefunctions with 
the same principal quantum number and decreases the 
overlap of states with different principal quantum num- 
bers. This expectation is reflected in the lattice sim- 
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Figure 1: The matrix elements for decay of T to r\b with the 
same principal quantum number, as a function of momen- 
tum. Shown are 0(v 4 ) (circles) and 0(v 6 ) (squares) action 
results with the leading 0(v 4 ) magnetic operator. Triangles 
include the relativistic correction to the magnetic operator in 
the 0(v 6 ) calculation. 
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Figure 2: The matrix elements for decay of an excited T 
to the rib ground state, as a function of momentum. Shown 
are 0(v 4 ) (circles) and 0(v 6 ) (squares) action results with 
the leading 0(v 4 ) magnetic operator. Triangles include the 
relativistic correction to the magnetic operator in the 0(v 6 ) 
calculation. 
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Table II: Three-point matrix elements from simultaneous fits to 12 correlation functions with different source-sink separation. 



T-t 3 


A (VP) 

A \\ 


a ( PV) 
A ^ ' 
A 1\ 


A ( PV) 

A * ' 

A m 


A (VP) 

A ■ ' 
A 21 


A (VP) 

A*' 
A 31 


A (VP) 

A ■ ' 
A 22 








momentum 








15 


0.916(2) 


-0.068(3) 


-0.050(3) 


0.071(5) 


0.070(7) 


0.863(67) 


19 


0.915(2) 


-0.068(2) 


-0.050(4) 


0.071(4) 


0.065(3) 


1.11(23) 


27 


0.916(2) 


-0.068(3) 


-0.051(6) 
momentum 1 


0.071(4) 


0.062(4) 


1.9(1.8) 


15 


0.907(2) 


-0.060(6) 


-0.047(5) 


0.082(5) 


0.072(5) 


0.799(68) 


19 


0.907(1) 


-0.062(6) 


-0.047(7) 


0.079(5) 


0.067(5) 


0.95(21) 


27 


0.907(2) 


-0.061(5) 


-0.048(8) 
momentum 2 


0.079(5) 


0.066(6) 


1.6(1.5) 


15 


0.877(1) 


-0.031(4) 


-0.040(5) 


0.103(6) 


0.083(6) 


0.791(60) 


19 


0.877(1) 


-0.031(4) 


-0.041(6) 


0.102(5) 


0.078(6) 


1.02(20) 


27 


0.878(2) 


-0.029(5) 


-0.039(7) 


0.100(5) 


0.068(6) 


1.0(1.6) 
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Figure 3: The matrix elements for decay of an excited Jj;, 
to the T ground state, as a function of momentum. Shown 
are 0(v 4 ) (circles) and 0(v 6 ) (squares) action results with 
the leading 0(v 4 ) magnetic operator. Triangles include the 
relativistic correction to the magnetic operator in the 0(v 6 ) 
calculation. 



ulations. Figure [T] compares transition amplitudes be- 
tween states with the same principal quantum number 
calculated with the 0(v 4 ) and 0(v 6 ) actions using only 
the leading magnetic 0(v A ) transition operator. The de- 
creased spin-dependent interaction at 0(v 6 ) shows up as 
an increase of the transition amplitude. The triangles 
in Fig. [1] show the effect of adding the relativistic cor- 
rection to the magnetic transition operator in the 0(v e ) 
calculation. The effect is very small for these transitions. 



Figure [2] shows the comparison of the calculations for 
the case of excited T decaying to the ground state r\h ■ 
Here going from 0(v 4 ) (circles) to 0(v 6 ) (squares) leads 
to a decrease of the amplitude in line with the notion 
that the overlap of wavefunctions is decreased. What is 
surprising is that including the relativistic correction to 
the magnetic spin-flip (triangles) leads to such a large ad- 
ditional decrease of the amplitude. For the decay of an 
excited r]b to the ground state of T the situation is more 
complicated. The effect of spin-dependent interactions 
on the wavefunctions leads to a relative negative sign 
of the transition amplitude compared to excited T de- 
cay. Figure [3] shows that changing the action from 0(v A ) 
(circles) to 0(v 6 ) (squares) results in a decrease in the 
magnitude of the amplitude but including the relativis- 
tic correction to the operator (triangles) has the opposite 
effect. 

Up to now we have discussed only the transition due to 
magnetic spin-flip operators. There are additional con- 
tributions that have to be considered, that is, due to 
the electric terms. These are the electromagnetic coun- 
terparts of the spin-orbit terms in @ and ([3]) (see also 
Eq. (2) in Ref. [8]). These terms are suppressed by a 
factor of (photon energy/Mo) relative to the magnetic 
terms at the same order in v 2 . They make no contri- 
bution to the transition matrix element at zero photon 
momentum. For vector to pseudoscalar transitions be- 
tween states with the same principal quantum number, 
where the mass difference is small, their effect is not im- 
portant. For excited state to ground state transitions 
they make a non-negligible contribution. 

Table [TTTI gives our final results for excited state decay 
amplitudes interpolated (for 2S decays) or extrapolated 
(for 3S decays) to the physical momentum. The values in 
the first two rows were obtained using an NRQCD action 
with terms up to the indicated order but only the lead- 
ing magnetic spin-flip transition operator. For the third 
and fourth rows all terms in the transition operator up to 
the order of those included in the NRQCD action were 
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Table III: Three-point matrix elements for excited states decays at the physical momentum. 





T(25) rj b [lS) 
611 MeV/c 


T(3S) -» r?6(lS) 
923 MeV/c 


r?t,(2S) T(15) 
524 MeV/c 


7M3S) -» T(1S) 
830 MeV/c 




0.088(5) 


0.084(7) 


-0.055(4) 


-0.041(5) 


0(v 6 ),a 


0.073(4) 


0.067(7) 


-0.033(3) 


-0.023(5) 


Complete 0(v 4 ) 


0.080(5) 


0.077(7) 


-0.050(4) 


-0.032(5) 


Complete 0(v 6 ) 


0.032(5) 


0.013(7) 


-0.059(4) 


-0.041(5) 


Experiment 


0.034(7) 


0.018(3) 







evaluated. The difference between the first and third 
rows is due to the 0(v 4 ) electric term. The difference 
between the second and fourth rows is dominated by the 
correction to the magnetic spin-flip coupling. In the last 
row, the transition amplitude inferred from the experi- 
mentally observed T(2S) and T(3S) decays [3, [H[ is 
given. Our best estimate, complete 0(v 6 ), is consistent 
with the experimental values. 

The excited state decays are so-called hindered tran- 
sitions and the amplitudes depend on the interplay of 
spin-dependent effects on the states, relativistic correc- 
tions and recoil (momentum) effects. The systematics of 
excited T and excited r\\> decays are different and exper- 
imental information on r/i, decays would provide a very 
stringent test of calculations. Now that there are exper- 
imental glimpses [HI, [l3| of r]b{2S) one may hope that a 
measurement of its radiative decay may be feasible in the 
not too distant future. 

In this brief report we have updated a previous cal- 
culation [l| of bottomonium radiative decay amplitudes 
by extending the lattice NRQCD action to include spin- 
dependent terms of 0(v 6 ) and by including contributions 
to the transition operators that are consistent with the 
action. These changes bring the calculated decay ampli- 



tudes into agreement with values determined from ob- 
served excited T decays. The NRQCD action and tran- 
sition operators were constructed using tree-level coef- 
ficients (tadpole improved for the action). Radiative 
corrections to the coefficients were not considered here. 
These have been calculated for some cases and notice- 
able effects were found [1, [l6j]. A further refinement of 
the calculation would have to apply such corrections in a 
complete and consistent way. 

Finally, we note that the large changes in the excited 
state decay amplitudes found in going from 0(v 4 ) to 
0(v e ) NRQCD may suggest that it would be benefi- 
cial to avoid nonrelativistic approximations altogether. 
The construction and use of relativistic lattice actions 
for heavy quarks is the subject of ongoing investigations 
[13, EH ■ The comparison of nonrelativistic and relativis- 
tic approaches to excited bottomonium radiative decays 
would be a worthwhile direction for future work. 
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gineering Research Council of Canada. 
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